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A direct incident fluence measurement method based on amorphous silicon electronic portal imaging device 
(a-Si EPID) has been developed for pretreatment verification of intensity-modulated radiation therapy (IMRT). 
The EPID-based incident fluence conversion method deconvolves EPID images to the primary response dis- 
tribution based on measured lateral scatter kernels in the EPID detector using Conjugate Gradient algorithm. 
The primary response is converted to the incident fluence based on measured fluence conversion matrix which 
corrects for off-axis position dependence of the a-Si EPID response and the “horn” beam profile caused by 
flatting filter. To verify feasibility and accuracy of this method, square fields of various sizes and two IMRT 
plans were delivered. The dose distributions computed based on EPID-derived incident fluence were compared 
with the measurement data. For all square field sizes except the smallest field (2 cm), the mean dose differences 
in cross-line dose profiles were within 1% excluding the penumbra region, and gamma passing percentages 
with a 2%/2mm criterion were about 99%. For two IMRT plans, the least gamma passing percentage for all 
eight IMRT fields was 98.14% with 2%/3 mm criteria. It can be concluded that our direct EPID-based incident 
fluence conversion method is accurate and capable of being applied to pretreatment dose verification in clinical 


routines. 
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I. INTRODUCTION 


Intensity modulated radiotherapy (IMRT) is a powerful 
technique to achieve a better dose conformity to the tumor 
volume and reduce the dose to normal tissues as much as 
possible. This is only achieved if the planned fluence is de- 
livered accurately at the treatment unit. However, under-dose 
of the tumor or an over-dose of the organs at risk would be 
caused by the following factors [1]: (a) the calculated leaf 
sequence does not result in the fluence pattern used by the 
treatment plan system for dose calculation; (b) the treatment 
plan is not correctly transferred to the accelerator; and (c) the 
treatment equipment is not functioning correctly. Pretreat- 
ment verification of the fluence delivered is essential to en- 
sure the accuracy of IMRT treatment. Nevertheless, no detec- 
tor can be used so far to directly measure the intensity fluence 
for radiotherapy. 

Electronic portal imaging device (EPID), an X-ray imag- 
ing detector attached to the accelerator, has gone on to replace 
traditional dosimetry (such as ion chamber, film, etc.) as pre- 
treatment dose verification tool, due to the advantages of fast 
image acquisition, high resolution, digital format and good 
dosimetric characteristics [2-5]. One way to use the EPID 
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for pretreatment verification is to compare measured image 
(raw or converted to dose-to-water) with predicted EPID im- 
age [6-8] or dose-to-water at the level of the EPID [9, 10]. 
But its limitation is that it cannot directly interpret the dis- 
crepancy between the planned doses with the delivered dose 
inside patient. A more beneficial way to use EPID is to com- 
pare 3D dose reconstructed with plan calculated inside pa- 
tient, provided that the EPID image can be converted to the 
delivered energy fluence used as input for dose reconstruc- 
tion. Warkentin et al. [11] deconvoluted the EPID raw image 
to 2D primary fluence distribution using scatter kernels gener- 
ated by Monte Carlo model of EPID. Van Elmpt et al. [12-14] 
extracted the energy fluence from a measured dose-to-water 
with EPID by a back-projection model. Zijtveld et al. [15] 
used an iterative algorithm to obtain the incident fluence from 
portal dose image. However, modelling the EPID by Monte 
Carlo accurately is difficult, as detailed structure of EPID 
is usually not available. In Refs. [12—15], for obtaining the 
fluence, the EPID images should be converted to portal dose 
distribution by EPID dose calibration model, a complicate job 
with easy dosimetry errors. 


In this paper, a method is proposed to directly convert the 
EPID image to the incident fluence. It has an effective correc- 
tion procedure for the field size and off-axis position depen- 
dence of the EPID response, based on measured scatter ker- 
nels and fluence conversion matrix, which are derived from a 
series of measurements. EPID grey images are deconvolved 
to the primary response distributions based on scatter kernels 
using Conjugate Gradient algorithm. The primary responses 
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are converted to the incident fluence. To validate accuracy 
of EPID-derived incident fluence, the dose computed based 
on the measured fluence is compared to measured dose for a 
large range of conditions. 


Il. MATERIALS AND METHODS 
A. Electronic portal imaging device (EPID) 


The EPID is of a PerkinElmer 1640L amporphous silicon 
type EPID (PerkinElmer Optoelectronics, Wiesbaden, Ger- 
many) mounted on a Elekta Axesse linac (Elekta AB, Stock- 
holm, Sweden). On top of the sensitive layer of the EPID, a 
3mm copper slab was mounted to remove low energy pho- 
tons, which cause over-response of EPID [13, 16]. 

All images were processed to correct for pixel-to-pixel sen- 
sitivity variations and dark current. In addition, image lag and 
ghosting artifacts affecting amorphous silicon EPIDs [4, 17], 
were corrected by using a ghosting correction factor G, simi- 
lar to the procedure described by Nijsten [18]. 


B. EPID-based incident fluence conversion model 
1. Determining the scatter kernels 


The response J (x, y) measured at a point (x, y) in the EPID 
consists of two parts: the primary response Ip(, y) of radia- 
tion reaching the EPID plane directly and the composite scat- 
ter response [s(a, y) from lateral scatter within the EPID and 
accelerator head. 


I(x, y) = Ip(z, y) + Is(z, y) (1) 


In our method, the scatter component is modeled as 
convolution of the primary response with a scatter kernel 
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(2) 
where, Ip(x’, y’) is the primary response at point (2’, y’), 
T(e—2',y-y’) = I(x — x’)? + (y — y')?]!? is the distance of 
point (x’, y’) to the measured point (x, y), K(r(z—2 y—y’), d) 
is related to the off-axis distance d (beam-softening) due to 
energy dependency of EPID response. 

In order to derive the scatter kernel, the primary response 
directly reaching to the EPID is calculated first. Since the 
scatter component Is (x, y) depends on field size, the primary 
component Ip(x, y) has to be field size independent. This fact 
can be used to estimate the primary contribution as follows. 
The total on-axis response of the EPID I(x, y) is experimen- 
tally determined as a function of field size f,, by irradiating 
the EPID with square fields of different sizes. For the limit of 
zero field size, the scatter component is negligible. Thus, the 
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primary response Ip(x, y) can be assessed by extrapolating 
the total response I(x, y) to a field size fs of 0cm x Ocm: 


Ip(x,y) © lim I(x, y). (3) 
fs—0 
From Eqs. (1) and (2), we have: 
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where the point (x = 0,y = 0) represent the center of the 
beam field, f(r(q,,)) is the off-axis primary response profile, 
and a; and b; are the kernel parameters. We introduced radial 
symmetric N Gaussian function to describe the scatter kernel. 
To accurately derive a; and b;, the scatter contributions for 
a set of square fields of different sizes are estimated by sub- 
tracting the primary response computed using Eq. (3). The 
scatter data represented as a function of field size are fitted 
by adjusting the kernel parameters. 


2. Primary response extraction 


Assuming that the incident intensity at the source-axis 
distance (SAD) of 100cm is composed of pencil beams of 
Fy x Fy finite size, the response for pencil beam (i, j) reach- 
ing directly to the EPID is the primary response Ip(i, j). We 
can derive that the total response value of the point (x, y) in 
the EPID is the sum of the primary response Jp(i, j) and the 
scatter component from other pencil beams: 
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where Ic(x,y) and Imay) are the calculated and measured 
response values at a point (x,y) of the EPID, respectively. 
Ez; x Ey is the size of EPID image used to reconstruct the 
primary response Jp(i,7), and K(x, y, {i,7}) is the scatter 
contribution for point (x,y) in the EPID from pencil beam 
(i, j). Equation (6) is a nonlinear equation with Fẹ x Fy vari- 
ances, so the process of solving it becomes an unconstrained 
optimization problem. The convergence criterion is based on 
the difference between the calculated and measured response 
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values on the EPID. An in-house developed Conjugate Gra- 
dient algorithm [19, 20] was used to fit all the pencil beam 
primary response Jp(i, 7) so as to minimize error o between 
the calculated and measured values for all the E, x Ey pixel 
points pooled together. This algorithm is efficient for solv- 
ing large-scale optimization problems and has been applied in 
ARTS [21, 22] (Accurate/Advanced Radiation Therapy Sys- 
tem) developed by the FDS Team, an interdisciplinary team 
devoting to the research and development of nuclear reactor 
physics [23], nuclear reactor materials [24-26], nuclear reac- 
tor engineering [27-29], numerical simulation and visualiza- 
tion [30], medical physics and environmental protection, etc. 


3. Incident fluence conversion 


In order to convert the primary response values Ip(i, 7) 
derived from acquired EPID image to the incident fluence 
(x,y), we have to determine the relationship between the 
primary response values of the EPID and the number of mon- 
itor units (MU), defined as a unit representing the linac output 
of the accelerator. 

In our work, the fluence conversion matrix covering the 
whole detector was established based on the measurement 
data, which accounted for a general off-axis variation in 
beam fluence profile and variation in detector response due 
to changes in off-axis energy spectrum of the beam. 

We used a CC13 ionization chamber with a 3.76 cm mini- 
phantom as build-up cap to measure the in-air off-axis ra- 
dio (OAR) profile at a source-axis distance (SAD) of 100 cm 
by scanning (IBA Dosimetry GmbH, Schwarzenbruck, Ger- 
many) the diagonals of a 25 x 25cm? reference field. OAR 
from each half of both diagonal beam profiles was used to 
calculate an average OAR(r(,,,.,)). Under the same field size, 
EPID images were acquired for the calibrated number of 
monitor units 100 MU. Then, the primary response value 
of the EPID for each pencil beam Jp,,,(x,y) could be re- 
constructed using Eq. (6). The fluence conversion matrix 
CU (a, y) is determined by dividing the primary response val- 
ues and the number of monitor units corrected by OAR. 


CU (2, y) = Ibai (2,Y)/[Gcai» OAR rey] (7) 


Finally, the incident fluence S(x, y) can be calculated for 
any field size and shape based on the primary response ex- 
tracted by Step B and the fluence conversion matrix estab- 
lished by Eq. (7). 


C. Experimental validation of EPID-derived incident fluence 


Since direct measurement of incident fluence is difficult, it 
is impossible to estimate the accuracy of EPID-derived inci- 
dent fluence by comparing EPID-derived fluence with mea- 
sured fluence. Instead, comparison of dose calculated us- 
ing EPID-derived incident fluence to measured dose was per- 
formed to verify the EPID-derived incident fluence. An ac- 
curate in-house developed dose calculation engine based on 
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Monte Carlo Finite-Size Pencil Beam Model [31, 32] was 
adopted for dose calculation. 


1. Open field verification 


For verifying the sensitivity and accuracy of the EPID- 
based fluence conversion procedure for different field sizes 
and numbers of monitor unit (MU), square field sizes of 
2cm x 2cm, 5cm x 5cm, 10cm x 10cm, 15cm x 15cm 
and 20cm x 20cm (100 MU) and monitor units of 2, 5, 
10, 50, 100 and 150 MU (10cm x 10cm) were measured. 
The absolute doses at the field center and the cross-line 
dose profiles were measured at a depth of 5cm using the 
CC13 chamber in the Blue Phantom water phantom system 
(IBA Dosimetry GmbH, Schwarzenbruck, Germany) with a 
source-to-surface distance (SSD) of 100 cm. 


2. Clinical cases verification 


For verifying feasibility of the incident fluence measure- 
ment procedure for clinical cases, two step-and-shoot IMRT 
treatment plans were performed. The IMRT-1 treatment plan 
was a four-field lung plan with a total of 96 segments. The 
IMRT-2 treatment plan was a four-field head and neck plan 
with a total of 80 segments. Two-dimensional (2D) dose dis- 
tribution was measured at 5cm depth in a solid water phan- 
tom using 2D-ARRAY seven29 dosimetry (PTW-Freiburg, 
Germany) with an SSD of 95cm. The 2D dose distribution 
for each field was evaluated by comparing with the measured 
dose distribution. 


Ill. RESULT AND DISCUSSION 
A. Scatter kernels 


Scatter kernels for off-axis positions of 0, 2, 4, 6 and 8cm 
were derived from measurements. For the determination of 
these scatter kernels, square fields of 2, 5, 10, 15 and 20cm 
widths were measured with EPID for each off-axis position. 
The primary responses for different off-axis positions were 
calculated using Eq. (3). The coefficients a; and b; of the 
scatter kernel were determined by a fitting procedure using a 
non-linear least squares method implemented in Matlab (Mat- 
lab 7.6.0, The Mathworks, Natick, USA) based on Eqs. (4) 
and (5). The fitting procedure demonstrated that the best ker- 
nel fits for each off-axis position could be obtained if three 
Gaussian distributions (VV = 3) were used to describe the 
kernels. The scatter kernels for different off-axis position are 
shown in Fig. 1. The fitted parameters of the scatter kernels 
are given in Table 1. 

Figure 2 shows that the Gaussian distributions of the scat- 
ter kernel for different off-axis positions are close to each 
other. The small differences are possibly due to the impact of 
changes in the off-axis energy spectrum and the photons inci- 
dent direction. The same EPID image was deconvolved to the 
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Fig. 1. Scatter kernels for off-axis position of 0-8 cm. 


TABLE 1. Fitted parameters of scatter kernels for different off-axis 
positions 


Off-axis 

distance ai a2 a3 bı bo b3 
(cm) (cm7?) (cm~?) (cm~?) (cm) (cm) (cm) 
0 0.02197 0.003422 0.4900 1.601 4.760 —0.2509 
2 0.02184 0.003370 0.4886 1.588 4.885 —0.2467 
4 0.01956 0.003015 0.4887 1.665 5.272 —0.2467 
6 0.02208 0.003514 0.4896 1.591 4.622 —0.2503 
8 0.02352 0.004250 0.4899 1.601 4.760 —0.2502 


primary responses using scatter kernels for different off-axis 
positions in turn. Then, the primary responses were compared 
and the differences were less than 1%. Thus, the assumption 
that the field size-dependent scatter kernel is invariant with 
off-axis position is reasonable, as long as no object is in the 
beam. 


B. Experimental validation of EPID-derived incident fluence 
1. Open field verification 


The measured and calculated absolute dose profiles in 
cross-line direction for monitor units of 2, 5, 10, 50, 100 
and 150 MU (10cm x 10cm) and field sizes of 2cm x 2cm, 
5cm x 5cm, 10cm x 10cm, 15cm x 15cm, 20cm x 20cm 
(100 MU) are shown in Fig. 2. Gamma values are analyzed 
over the area within 5% isodose line of maximum dose using 
the criteria of 2%/2 mm. 

The mean dose differences of 2 and 5 MU are slightly 
larger than those of monitor units of > 10 MU, but they were 
still within 1.5%, excluding the penumbra region. For the 
gamma values, 95.5% of 2 MU and 98.7% of 5 MU were less 
than 1, with the mean gamma values being 0.627 and 0.569, 
respectively. For the other monitor units of > 10 MU, all 
gamma passing percentages were more than 99% and corre- 
sponding mean gamma values were less than 0.4. The larger 
dose differences for small numbers of monitor units can be 
possibly explained by the instable output of the linac. 

Except for 2cm x 2cm field, the mean dose differences 
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Fig. 2. The absolute dose profiles for different number of monitor 
units in 10cm x 10cm (a) and different field sizes in 100 MU (b). 
Solid lines, measured dose; dotted lines, computed dose based on 
the EPID-derived fluence. 


for all field sizes were less than 1% (excluding the penumbra 
region); gamma passing percentages were about 99% with 
mean gamma values less than 0.4. The dose error at the center 
of the field for the 2cm x 2cm field size was 2.7% and the 
corresponding gamma passing percentage was 69.2%. The 
larger dose error for 2 cm x 2 cm field size may be caused by 
a larger uncertainty due to the small field size relative to the 
ion chamber volume. 

The mean dose differences become slightly larger if the 
penumbral regions are included. But the gamma analysis 
indicates that the calculated doses based on EPID-derived 
fluence and the measured doses agree quite well within 5% 
isodose line of maximum dose, including the penumbral re- 
gions. This implies that the calculated dose profiles have a 
steeper dose falloff in the penumbra than the measured dose 
profiles due to blurring of the measurement by the ionization 
chamber. 


2. Clinical cases verification 


Two clinical step-and-shoot IMRT plan, consisting of eight 
fields totally, were measured with EPID and 2D ionization 
chamber array, respectively. VeriSoft5.1 software (PTW) was 
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Fig. 3. (Color online) Isodose lines of the measured dose (solid lines) and computed dose (dashed lines) based on the EPID-derived fluence 
for pretreatment verification of an IMRT field (a) and histograms of the gamma distribution (b). 


used to compare the measured dose distributions with the cal- 
culated dose distribution based on EPID-derived fluence. Fig- 
ure 3 shows the results for one field consisting of 26 segments 
as an example. The isodose lines (Fig. 3(a)) demonstrate very 
good agreement. Figure 3(b) shows the gamma histogram 
for this field. The gamma distributions with the criteria of 
2%/3cm for all eight IMRT fields were calculated over the 
area within the 5% isodose line of the maximum dose, which 
disregard the regions outside the field. The gamma values for 
all eight IMRT fields are given in Table 2. 


TABLE 2. Analysis of gamma-values for all IMRT fields using 
VeriSoft software (PTW) 


Plan Field Ymean Ymax Py<ı( %) 
IMRT1 1 0.385 1.094 99.09 

2 0.388 1.112 99.14 

3 0.408 1.102 98.78 

4 0.396 0.967 100.00 
IMRT2 1 0.363 0.952 100.00 

2 0.415 1.095 98.48 

3 0.396 0.814 100.00 

4 0.425 1.132 98.14 


All gamma distributions had a mean value below 0.43. The 
least gamma passing percentage (Py<1) for all IMRT fields 
was 98.14%. The gamma passing percentages for three fields 
reached 100%, implying that all the measured points satisfied 
the chosen gamma criteria. Excellent agreement between 
computed dose based on EPID-derived incident fluence and 
measured dose indicates that the incident fluence obtained by 
our method is accurate. 

The reconstructed incident fluence can be used as input 
to TPS, so as to generate a 3D dose distribution using the 
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patient’s CT data. It can provide the cumulative errors in 
a 3D dose distribution in the patient from all beams in the 
IMRT plan to help the clinician to assess the potential clinical 
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IV. CONCLUSION 


In this article, we have developed an accurate EPID-based 
incident fluence conversion procedure to directly convert 
EPID images to the incident fluence. For this purpose, a 
conversion model has been investigated that corrects for the 
field size and the off-axis position dependence of the EPID 
response based on measured scatter kernels and fluence con- 
version matrix, without the need of input of beam characteris- 
tics and detailed analysis of EPID construction. EPID-derived 
incident fluences were evaluated for different fields such as 
open fields and IMRT fields to determine the accuracy and 
applicability of the EPID-based incident fluence conversion 
model. Gamma evaluation criteria were satisfied for all mea- 
surements. Using our EPID-based incident fluence conver- 
sion method, highly accurate incident fluence can be obtained 
which are promising for pretreatment 3D dose verification. 
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